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SECT  I  ON  1 


INTRODUCTION 

This  report  documents  modifications  and  additions  incorporated  into  the 
ABRES  Shape  Change  Code  (ASCC80)!  to  create  a  planar  two-dimensional  version 
of  the  axisymmetric  computer  code.  This  planar  code  predicts  convective  heat 
transfer  and  in-depth  conduction  for  swept  wings  in  supersonic  flow.  These 
planar  modifications  to  ASCC80  were  developed  under  the  Aerodynamic  Heating 
Computations  for  Projectiles  program.  The  overall  objectives  for  the  program 
were  threefold: 

1.  Modify  the  in-depth  heat  conduction  pacVage  to  improve  ASCC's 
capabilities  to  handle  slender  mul timaterial  configurations 

2.  Extend  the  developments  of  planar  ASCC  modifications  to  predict 
heating  of  swept  fin  configurations  to  include:  (a)  turbulent  flow 
on  swept  wings;  (b)  2-D  shock  shape;  and  (c)  improved  in-depth  heat 
conduction  routines 

3.  Develop  an  interactive  computational  grid  develping  routine  to 
simplify  the  procedure  for  specifying  body  configurations  and 
developing  computational  grids  for  ASCC 

The  modifications  made  to  ASCC80  covering  the  second  objective  are 
documented  in  Volume  II  of  this  report.  Volumes  I  and  III  of  this  report 
document  the  work  related  to  Objectives  1  and  3,  respectively.  In  this 
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document,  the  updated  ASC  code  Is  referred  to  as  PINRASCC,  and  the  updated 
code  associated  with  Objective  1  Is  referred  to  as  BRLASCC. 

Technical  discussion  of  the  PLNRASCC  modifications  Is  presented  In 
Section  2.  Section  3  Is  devoted  to  a  discussion  of  input  and  output. 
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SECTION  2 


TECHNICAL  DISCUSSION 

2.1  VISCOUS  AND  INVISCID  FLOWFIELD  MODIFICATIONS 

Viscous  and  inviscid  flow  models  for  the  ASCC80  version  of  the  ABRES 
Shape  Change  Code^  have  been  modified  to  solve  a  planar  two-dimensional  model 
for  flow  over  swept  wings.  The  theoretical  basis  for  these  swept  wing 
modifications  is  described  in  a  report  by  Suchsland.^  Implementation  into  the 
PLNRASCC  computer  program  is  taken  directly  from  that  report.  In  addition, 
two  new  capabilities  have  been  added  to  the  fluid  flow  modeling  in  the 
computer  code.  The  first  is  the  ability  to  model  two-dimensional  shock  shapes 
and  include  these  effects  In  the  boundary  layer  calculation.  Suchsland's 
version  of  the  code  could  only  model  shock  shapes  for  axisymmetric 
configurations.  The  second  added  capability  is  the  use  of  curve  fitted 
pressure  correlations  for  planar  geometries.  The  earlier  version  of  the 
inviscid  flow  model  was  restricted  to  axisymmetric  body  shapes. 

The  planar  shock  shape  modifications  to  the  code  are  a  series  of  curve 
fits  to  computed  results  from  a  two-dimensional  version  of  the  RAZZ  I 
computer  code.  The  RAZ2IB  code  is  a  general  inviscid  flow  solver  for 
highspeed  flight  configurations.  Pressure  distributions  were  computed  for  a 
series  of  cylinder-wedge  two-dimensional  wings.  Flow  conditions  for  these 
test  cases  spanned  a  range  of  Mach  numbers  (f’U  =  1.75,  2.0,  3.0,  4.0,  and 
6.0),  and  also  a  range  of  aft  wedge  angles  (ew  =  0.0°,  2.5°,  and  5,0°).  The 
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Figure  2-1.  Sketch  of  swept  wing  geometry 


computed  results  for  shock  angle,  0W,  versus  X/R,  were  fitted  with  a  least 
squares  sixth  degree  polynomial  for  each  combination  of  wedge  angle  and  Mach 
number  and  implemented  into  the  BRLASC  code.  Internally,  these  curve  fits 
are  converted  into  tabular  form.  Linear  interpol ation  in  both  wedge  angle  and 
Mach  number  is  used  to  create  the  shock  shape  table  for  each  new  swept  wing 
geometry.  Note  that  in  the  case  of  a  swept  wing,  the  required  shock  shape  is 
computed  using  the  component  of  freestream  Mach  number  that  is  normal  to  the 
leading  edge  of  the  wing. 

The  planar  invlscid  pressure  correlations  are  also  implemented  in  the 
form  of  curve  fits  to  inviscid  flow  calculations.  Results  from  the  RAZZ  I B 
code  were  obtained  for  the  same  ranges  of  Mach  numbers  and  wedge  angles 
described  above.  Curve  fits  for  P/P0  versus  X/R-j  were  produced  using  a 
combination  of  sixth  degree  polynomial  and  exponential  least  squares 
functions.  Linear  Interpolation  is  used  in  both  wedge  angle  and  Mach  number 
in  order  to  produce  a  pressure  ratio  for  each  boundary  layer  integration 
point.  As  was  noted  above  for  the  shock  shape  predictions,  it  is  the  swept 
wing's  normal  component  of  Mach  number  that  determines  where  to  interpolate  in 
the  pressure  curve  fits. 


The  new  planar  pressure  correlations  are  only  used  on  the  aft  wedge 
portion  of  the  wing.  The  existing  ASCC80  pressure  correlations  use  modified 
Newtonian  theory  to  predict  the  pressure  distribution  on  the  nosetip  of  an 
ax  i  symmetric  configuration.^  Newtonian  theory  applies  to  planar  geometries  as 
well  as  axisymmetric  ones.  Thus  the  original  pressure  correlations  were  not 
altered  in  the  nosetip  region. 

Planar  input  modifications  to  the  BRLASC  code  are  described  in 
Section  3.  This  section  should  be  used  in  addition  to  the  user's  manual  in 
Volume  I  of  this  report  in  order  to  run  the  new  code  for  planar  swept  wings. 
These  Input  modifications  enable  PLNRASCC  to  be  considerably  more  versatile 
and  easier  to  use  than  Suchsland's  original  version  of  the  code. 

2.1.1  Resul ts 

This  section  gives  convective  heat  transfer  results  from  PLNRASCC  and 
compares  them  with  experimental  data  from  a  number  of  different  swept  wing 
configurations.  These  experiments  were  chosen  to  test  the  code  for  its 
ability  to  compute  heat  transfer  in  both  laminar  and  turbulent  boundary 
layers . 

Figures  2-2  through  2-5  show  predictions  of  the  swept  wing  data  of 
Stainback.5  This  experiment  consisted  of  a  60°  swept  delta  wing  in  a 
supersonic  flow  at  several  different  tunnel  stagnation  pressures.  The 
freestream  Mach  number  for  these  data  is  4.95.  Laminar  flow  conditions  exist 
for  the  entire  run  length  of  the  wing.  The  "new  prediction"  in  Figures  2-2 
through  2-5  was  produced  using  PLNRASCC.  The  "old  prediction"  is  taken  from 
Suchsland.2  The  differences  between  these  two  codes  lie  in  the  formulation  of 
shock  shape  and  inviscid  pressure  results  for  planar  geometries.  In 
Suchsland's  prediction,  he  assumes  a  normal  shock  at  the  leading  edge  of  the 
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Figure  2-2.  Heat  transfer  predictions  of  the  experimental  data  of  Stainback- 
Rj  =  0.25  inch,  P0  =  428  psig,  T0  =  460°F 


NORMALIZED  RUN  LENGTH  (s/r) 

Figure  2-3.  Heat  transfer  predictions  of  the  experimental  data  of  StainbackS 
Ri  =  0.25  inch,  P0  =  223  psig,  T0  =  44}°F 


Figure  2-4.  Heat  transfer  predictions  of  the  experimental  data  of  Stainback^ 
Rj  «  0.25  inch,  P0  =  109  psig,  T0  =  432°F 
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Figure  2-5.  Heat  transfer  predictions  of  the  experimental  data  of  Stainback^ 
R,  =  0.25  inch,  P0  *  65  psig,  T0  =  437°F 
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wing.  Also,  the  pressure  ratio  ( P /P0 ) »  obtained  from  the  ax 'symmetric 
ASCC80  correlation. 

The  results  of  Figures  2-2  through  2-5  indicate  that  PLNRASCC  and 
Suchsl and ' s  version  of  ASCC  produce  very  similar  results  for  these  cases. 

There  are  two  reasons  for  this.  First,  the  normal  component  of  Mach  number  to 
the  wing  leading  edge  is  relatively  low,  (2.47).  This  means  that  the 
influence  of  the  shock  shape  on  the  boundary  layer  will  also  he  small.  For 
higher  Mach  numbers,  the  shock  shape  will  have  a  more  significant  effect  on 
entropy  swallowing  in  the  boundary  layer.  The  second  reason  for  the 
similarity  in  predictions  is  that,  for  zero  wedge  angle,  the  pressure 
predictions  for  the  axisymmetric  and  two-dimensional  correlations  will  be 
similar.  This  is  not  necessarily  the  case  for  higher  wedge  ang'es. 

Figures  2-6  through  2-9  show  predicted  and  experimental  -esults  for  the 
data  of  Murray  and  Stallings.6  Their  experiment  tested  both  60"  and  ?UU  swept 
wings  in  a  w’nd  tunnel  with  a  range  of  freestream  Mach  numbers  and  tunnel 
stagnation  pressures.  The  aft  wedge  angle  was  zero  for  all  of  these  cases.  A 
boundary  layer  trip  was  used  at  S  =  0.637  cm  from  the  wing  leading  edge  in 
order  to  provide  for  turbulent  flow  over  the  wing. 

For  the  70°  swept  wings  tested  by  Murray  and  Stallings,  the  normal 
component  of  Mach  number  to  the  wing  leading  edge  is  too  low  (<  1 . 75)  to  be 
covered  by  ttie  pressure  and  shock  shape  correlations  described  in  Section  2.1 
of  this  report.  Therefore,  predicted  results  are  only  presented  for  cases 
with  a  60°  sweep  angle.  These  predictions  were  made  with  PLNRASCC.  The 
calculation  was  done  with  specified  transition  to  turbulent  flow  at  the 
boundary  layer  trip.  Transitional  heating  was  modeled  in  this  calculation  by 
specifying  NREYCR  =  +4  in  Input  Table  1  of  the  input  data. 
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Predicted  results  cannot  be  expected  to  agree  with  the  data  near  the 
boundary  layer  trip  because  the  boundary  layer  requires  a  finite  length  after 
the  trip  until  it  becomes  fully  turbulent.  This  boundary  layer  recovery 
distance  Is  clearly  seen  in  the  experimental  data  of  Figures  2-6  through  2-9. 
The  computer  model,  on  the  other  hand,  can  change  abruptly  from  laminar  to 
turbulent  flow. 

Agreement  between  the  PLNRASCC  prediction  and  the  It  =  3.71 
experimental  data  in  Figures  2-6  and  2-7  is  similar  to  what  is  found  by  Murray 
and  Stallings.*’  They  present  predicted  results  that  were  obtained  from  the 
"strip  theory"  method  of  Van  Driest. ^  For  the  higher  Mach  number  cases  in 
Figures  2-7  and  2-8,  PLNRASCC  underpredicts  the  heat  transfer  on  the  aft 
portion  of  the  wing. 

Murray  and  Stallings®  estimate  the  uncertainty  in  their  Measured  heat 
transfer  coefficients  to  be  10  percent  for  h  >  306  15  percent  for 

20  J/m^-s-K  <  h  <  306  J/m^-s-K,  and  20  percent  for  h  <  20  J/m^-s-K.  The 
predicted  heat  transfer  results  from  PLNRASCC  are  close  to  being  within  these 
uncertainty  ramies.  Changes  introduced  to  the-  experimental  boon  iary  layer  by 
the  presence  of  the  trip  may  account  for  the  remainder  of  the  di screpancies. 

It  should  be  noted  that  in  all  cases,  the  basic  shapes  of  the  PLNRASCC  heat 
transfer  predictions  show  good  agreement  with  the  experimental  dita. 

A  final  ‘omparison  between  PLNRASCC  am  experinvntal  da f -  Is  shown  In 
Figures  2-10  through  2-12.  These  experiments  were  conducted  by  hunt  et  al  .fi 
They  used  temperature-sensi tlve  paint  on  a  60°  swept  wing  In  order  to  measure 
the  heat  transfer  rates.  High  uncertainties  are  usually  associated  with  this 
phase  change  paint  technique,  although  no  specific  values  are  provided  by  Hunt 
et  al .8  The  aft  wedge  angle  was  zero  for  all  of  their  experiments.  These 
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Figure  2-10.  Heat  transfer  predictions  of  the  experimental  data  of  Hunt 

©  et  al  .«  Ri  =  0.5  inch,  =  7.81 , 

£>  Re  =  0.92  x  10$  (based  on  leading  edge  diameter),  a  =  60° 
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Figure  2-11.  Heat  transfer  predictions  of  the  experimental  data  of  Hunt 
et  al .3  R1  =  0.5  Inch,  K,  *  7.94, 

Re  =  2.6  x  10$  (based  on  leading  edge  diameter),  A  =  60° 


Figure  2-12.  Heat  transfer  predictions  of  the  experimental  data  of 
Hunt  et  al .8  Rj  =  0.5  Inch,  =  7.98, 

Re  =  9.3  x  10^  (based  on  leading  edge  diameter),  a  =  60° 

experiments  had  freestream  Mach  numbers  close  to  8.0  and  a  range  of  Reynolds 
numbers. 

Results  from  the  low  Reynolds  number  case  in  Figure  2-10  are  in  fairly 
good  agreement  with  the  experimental  data.  The  largest  discrepancies  occur 
near  the  leading  edge  of  the  wing,  laminar  flow  was  observed  experimentally 
over  the  entire  wing,  and  the  PLNRASCC  result  represents  a  laminar  boundary 
layer  calculation. 

Figure  2-11  gives  heat  transfer  results  for  a  higher  Reynolds  number 
cjse.  Again,  the  PLNRASCC  prediction  Is  In  fairly  good  agreement  with  the 
data.  Laminar  flow  was  observed  experimentally  over  the  full  length  of  the 
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wing,  so  a  laminar  boundary  layer  calculation  was  performed  with  the  PLNRASC 
code. 

Heat  transfer  results  for  the  highest  Reynolds  number  case  of  Hunt 
et  al.®  are  shown  In  Figure  2-12.  Turbulent  flow  was  observed  on  the  surface 
of  the  wing,  but  it  is  difficult  to  determine  exactly  where  boundary  layer 
transition  takes  place.  Hunt  et  al.®  concluded  from  their  data  that  the  flow 
was  turbulent  on  the  cylindrical  nose  of  the  wing.  It  apparently  laminarized 
as  It  expanded  around  the  nose  and  onto  the  flat  plate.  This  laminarized 
boundary  layer  then  went  through  a  transition  to  turbulent  flow  at 
approx imately  S  =  1.2  Inches  on  the  flat  plate  portion  of  the  wing. 

The  PLNRASC  code  has  no  mechanism  for  modeling  this  sort  of  a 
1  arninarizing  boundary  layer.  An  attempt  was  made,  however,  to  test  the 
PLNRASC  code  against  this  experimental  data.  Two  separate  computer 
predictions  were  made.  The  first  was  a  fully  laminar  calculation.  The  second 
prediction  used  a  fixed  transition  location  at  S  =  0.2  inch  from  the  leading 
edge  of  the  wing.  Figure  2-12  Indicates  that  the  turbulent  prediction  shows 
better  agreement  with  the  experimental  data,  particularly  In  the  region  away 
from  the  nose  of  the  wing.  Neither  prediction  does  well  at  the  leading  edge 
of  the  wing. 

It  should  be  noted  that  a  fully  turbulent  PLNRASCC  calculation  was  also 
attempted  for  this  case.  Abrupt  transition  to  turbulence  was  specified 
immediately  after  the  laminar  series  solution  that  starts  off  the  calculation. 
This  caused  the  boundary  layer  properties  Influence  coefficient  on  Stanton 
number  to  change  abruptly.  The  effect  of  this  coefficient  Is  very  Important 
at  the  leading  edge  of  the  wing.  The  result  was  that  the  turbulent  heat 
transfer  rate  dropped  below  the  corresponding  laminar  heat  transfer  rate  for 
the  first  few  integration  points.  This  was  clearly  an  unrealistic  result,  so 


20 


the  transition  location  was  moved  farther  back  on  the  nose  until  this  problem 
was  no  longer  encountered. 

Hunt  et  al.®  also  described  a  calculation  method  that  they  used  to 
predict  the  heat  transfer  rates  In  their  experiments.  This  finite  difference 
calculation  method  included  a  spanwlse  momentum  equation  for  predicting  swept 
wing  cases  with  large  crossflows.  The  PLNRASCC  heat  transfer  predictions  in 
Figures  2-10  and  2-11  agree  quite  well  with  the  predictions  that  Hunt  et  al . 
present  In  their  paper.  For  the  highest  Reynolds  number  case  shown  in 
Figure  2-12,  the  laminar  PLNRASCC  prediction  is  in  good  agreement  with  a 
laminar  prediction  that  is  given  in  the  Hunt  et  al  .8  paper. 

2.1.2  Discussion  and  Conclusions  on  Flowfield  Modifications 

For  laminar  flow  at  relatively  low  supersonic  Mach  numbers,  PLNRASCC 
does  an  adequate  job  of  predicting  heat  transfer  on  swept  wings.  This  is 
evidenced  in  the  results  of  Suchsland2  and  also  in  the  results  presented  in 
Figures  2-2  through  2-6  of  this  report.  The  limitations  of  this  method  for 
laminar  flow  are  guided  by  the  basic  assumptions  that  were  used  in  the 
formulation  of  the  swept  wing  integral  equations  of  Suchsland.2 

The  PLNRASC  code  has  been  demonstrated  to  be  capable  of  predicting 
turbulent  boundary  layer  heat  transfer  over  swept  wings  in  supersonic  flow. 

The  predicted  results  shown  in  Figures  2-6  through  2-12  appear  to  be  quite 
reasonable  considering  the  uncertainties  associated  with  the  experimental 
data. 

One  area  that  needs  to  be  investigated  further  Is  the  modeling  of 
boundary  layer  transition  on  swept  wings.  Transition  criteria  based  on 
momentum  thickness  Reynolds  number  cannot  be  easily  applied  to  PLNRASCC  for 
swept  wing  cases  with  cross-flow.  This  Is  because  PLNRASCC  computes  the  flow 
in  a  two-dimensional  direction  that  is  normal  to  the  leading  edge  of  the  wing. 


Reynolds  numbers  are  formed  in  the  code  by  using  a  vector  component  of 
boundary  layer  edge  velocity.  Thus,  the  boundary  layer  momentum  thickness 
Reynolds  number  that  is  calculated  by  PLNRASCC  has  a  questionable  physical 
interpretation  regarding  boundary  layer  transition.  It  is  suggested  that 
users  of  PLNRASCC  supply  a  fixed  transition  location  to  the  code  that  is 
obtained  from  a  procedure  that  is  appropriate  to  three-dimensional  boundary 
1 ayers . 

2.2  PLANAR  IN-DEPTH  HEAT  CONDUCTION  MODIFICATIONS 

All  of  the  modifications  made  to  BRLASCC  have  been  incorporated  into 
PLNRASCC,  and  the  documentation  for  these  changes  can  be  found  in  Volume  I  of 
this  report.  However,  BRLASCC  is  an  axisymmetrlc  shape  change  code;  hence,  the 
heat  conduction  equation  formulations  used  In  PLNRASCC  must  be  modified  to 
correctly  model  the  planar  nature  of  the  problems  of  interest. 

2.2.1  Implicit  Grid  Modifications 

The  conduction  equation  in  the  moving  orth  ional  coordinate  system 
(implicit  grid)  under  the  axisymmetrlc  assumption  (a/ar  =  0)  Is: 
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(1) 


where 


Cp  =  specific  heat 

r0  =  body  circumferential  radius  of  curvature 
rb  =  r0  +  r»c  os  (9) 

rc  =  local  streamwlse  radius  of  curvature 
k  =  thermal  conductivity 
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p  ■  density 

n  «  surface  normal  recession  rate,  n  *  -r 
T  =  temperature 
t  =  time 

0  =  angle  between  normal  to  local  surface  and  axis  of  symmetry 
s  =  streamwlse  distance  along  body 

r  «  distance  normal  to  body  surface  at  s,  measured  from  surface 
Under  the  planar  assumption,  the  conduction  equation  takes  the  form: 
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The  finite-difference  equations  In  PLNRASCC  have  been  modified  to 
reflect  the  change  in  the  differential  equation,  setting  rb  to  unity 
wherever  the  body  radius  of  curvature  appears. 

2.2.2  Explicit  Grid  Modifications 

Making  use  of  the  axlsymmetrlc  nature  of  the  problems  of  Interest  the 
conduction  equation  in  BRLASCC  utilizes  a  fixed  cylindrical  coordinate 
system  (explicit  grid)  and  Is  given  by: 


pC 


3T 

p  TE 


(3) 


PLNRASCC  however  Is  not  axlsymmetrlc,  and  the  conduction  equation  takes 
the  form  of  the  two-dimensional  Cartesian  coordinate  heat  conduction  equation 
and  Is  given  by: 
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r  3T  .  3  /  3 T \  ,  3  /  3 T \ 

pCp  3t  ‘  3x  [K  3x  /  a.y  \*  3y ) 


(4) 


The  finite-difference  equations  in  PLNRASCC  have  been  modified  to 
reflect  the  change  in  the  differential  equation. 

2.2.3  Validation  of  Planar  Conduction  Modifications 

Analytical  solutions  presented  in  graphical  form  are  available  from 
Heisler^  for  the  transient  temperature  distribution  of  multidimensional 
systems.10  The  solution  for  the  transient  temperature  distribution  of  an 
infinite  copper  rectangle  4  inch  x  8  inch  was  used  to  check  the  planar 
conduction  modifications  made  in  PLNRASCC.  Figure  2-13  illustrates  tne 
problem. 

The  rectangle  is  at  a  uniform  initial  temperature  of  530°R.  The 
surface  temperature  is  suddenly  raised  to  800°R.  The  temperature  versus  time 
at  Point  A  was  calculated  using  the  Heisler  charts  and  compared  with  PLNRASCC. 
The  results  are  shown  in  Figure  2-14.  The  agreement  between  PLNRASCC  and  the 
Heisler  chart  solution  is  excellent.  Only  at  5  s  does  the  solution  from 
PLNRASCC  depart  significantly  from  the  Heisler  chart  solution.  This  is  due  to 
the  fact  teat  the  Fourier  modulus  at  5  s  is  less  than  0.2,  where  the  Heisler 
charts  become  invalid. 
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TEMPERATURE  (DEG  R) 

640.00  680.00  720.00  760  00  800.00 


BRLASCC  2-D  PLANAR  CONDUCTION 
UERSUS  HEISLER  CHART  SOLUTION 


A 


1 - 1 - 1 - 1 - 1 - 1 

0.00  10.0  20.0  30.0  40.0  50.0 

TIME  ( SEC  ) 


Figure  2-14.  Comparison  between  PLNRASCC  conduction  solution  and  Heisler 
chart  solution  for  an  Infinite  rectangle 
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INPUT  AND  OUTPUT 


3.1  INPUT  MODIFICATIONS  TO  THE  BRLASCC  COMPUTER  CODE 

Two  of  the  input  tablet  to  the  BRLASC  code  require  modification  in 
order  to  run  the  planar  version  of  the  code  for  swept  wings.  The  conventional 
forms  of  these  tables  are  described  in  Volume  I  of  this  report.  The 
modifications  that  are  described  below  should  be  implemented  in  conjunction 
with  this  reference. 

TABLE  7:  Surface  Data 

The  modification  to  this  table  consists  of  an  additional  surface  table 
that  describes  the  geometry  of  the  swept  wing  and  also  the  specified  option 
for  computing  the  inviscid  flowfleld.  Whenever  a  swept  wing  configuration  is 
to  be  computed,  the  following  subtable  must  be  included  in  the  surface  data. 
This  subtable  takes  the  following  form: 


Col umns 

Format 

Data 

Units 

1  -  5 

15 

Enter  4  (Swept  wing  configuration 
subtable) 

-- 

6-10 

15 

Option  for  pressure  calculation 
=  1  :  Use  2-D  cylinder/wedge 
correlations 

/  1  :  Do  not  use  2-D  cylinder/ 
wedge  correlations 

1  -  10 

F10.5 

Sweep  angle  of  the  wing  measured 
from  a  normal  to  the  flow  direction 

degree 

11  -  20 


F10.5 


Aft  wedge  angle 


degree 


k+2  1-5  15  Enter  -1  if  other  surface  tables 

follow;  +1  for  last  surface  table 

Note  that  if  the  2-0  cylinder/wedge  correlations  for  pressure  are  not  desired, 
the  user  has  a  number  of  options  for  specifying  pressure  that  are  described  in 
Volume  I  of  this  report.  These  include  two  ways  of  specifying  tabular  values 
for  pressure  as  well  as  using  the  axisymmetric  pressure  correlations  that  are 
built  into  BRLASCC  and  ASCC80. 

TABLE  8:  Shock  Shape  Data 

This  table  gives  the  user  the  option  of  whether  or  not  to  use  the  2-D 
shock  shape  correlations  for  swept  wing  configurations.  The  modification  for 
input  to  this  table  is  the  addition  of  another  shock  shape  flag.  This 
modification  is  described  below. 

Columns  Format  Data  Uni  ts 

1-2  12  Enter  08  (table  number) 

1-5  15  ISHFIG  --  shock  shape  flag 

1-  Shock  angle  given  as  function 
of  y  coordinate 

2-  Shock  angle  given  as  function 
of  dimensionless  y  coordinate 
(y/R  ),  where  R  is  the  nose 
radius. 

3-  Flag  to  use  the  2-D  planar  shock 
shape  correlations 

The  remainder  of  this  table  is  identical  to  the  one  in  Volume  I  of  this 
report.  Note  that  if  TABLE  8  is  not  specified  in  the  program  input  data,  the 
code  will  use  a  shock  shape  that  Is  generated  from  the  axisymmetric 
correlations  that  are  built  in  to  BRLASCC  and  ASCC80. 


Card  No. 
1 
2 


SAMPLE  PROBLEM  1 


Swept  Wing  Configuration 
Boundary  Layer  Solution 


3.2  SAMPLE  PROBLEM 


This  subsection  Illustrates  two  sample  cases  for  the  PLNRASC  code.  The 

5 

first  sample  case  Is  a  prediction  of  a  swept  wing  experiment  from  Stainback. 
Heat  transfer  results  from  this  case  are  presented  in  Figure  2-1. 

The  complete  input  data  for  this  case  is  presented  followed  by  selected 
portions  of  the  output  file.  The  user  should  note  that  the  boundary  layer 
calculation  proceeds  in  a  direction  that  Is  normal  to  the  leading  edge  of  the 
wing.  Output  quantities  such  as  boundary  layer  edge  velocity  represent  vector 
components  of  the  total  flowfiald  quantities. 

The  second  sample  problem  is  a  planar  calculation  of  the  sample  problem 
found  in  Volume  1  of  this  report.  The  axisymmetric  projectile  of  Volume  I  has 
been  made  into  a  planar  60°  swept  wing.  A  complete  listing  of  the  input  data 
is  presented  as  well  as  selected  output  data. 
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Sample  Problem  1  Output 


PLANAR  VERSION  BRL  IMPROVED  ABRES  SHAPE  CHANGE  CODE  (PLNARASCC) 


input 


BRL/ASCC  TWO-D I  MENS IONAL  TEST  CASE 
PLANAR  SWEPT  WING 


8/1 0 


—  general  program  r  lags  — 


(ENVIRONMENT  FlAG) 

LC  - 

2 

(SHAPE  CHANGE  F LAG) 

ISS  - 

2 

(OUTPUT  PRiNT  FLAG) 

IPRNT  - 

2 

(TRANSITION  CRITERIA  FLAG) 

NREYCR  - 

2 

(BOOY  ANGLE  OEFN  FLAG) 

IRON  - 

0 

(CARBON  TRANS  CRIT.  FLAG) 

ICARB  - 

1 

(NOSE  SHAPE  MODIFICATION  FLAG) 

IMOO  - 

0 

-  TIME  INCREMENT  INFORMATION  - 


INITIAL  TIME  (SEC)  0  0000  FINAL  TIME  (SEC)  0  0000 

OUTPUT  INTERVAL  -  1  0000  SEC  FROM  INITIAL  TIME  UNTIL  2  0000  SEC 

OUTPUT  INTERVAL  «  1.0000  SEC  FROM  2  0000  SEC  UNTIL  FINAL  TIME 


BOUNDARY  LAYER  ONLY  SOLUTION 

COMPUTATION  TIME  STEPS  SET  EQUAL  TO  SPECIFIED  ENVIRONMENT  TIME  STEPS 


-  WIND  TUNNEL  ENVIRONMENT  - 

FREESTREAM  MACH  NO  «=  4  95 

TIME 
(SEC) 

0.000 
1  .000 
2.000 
3.000 


TOTAL  PRE5SURE 
(PS1A) 

65  00 
109  00 
223  00 
420  00 


TOTAL  TEMPERATURE 
(DEG  E ) 
437.00 
432.00 
441  00 
460.00 


PLANAR  VERSION  SRL  IMPROVED  A PRES  SHAPE  CHANCE  COOE  (PLNARASCC) 
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planar  version  brl  improved  abres  shape  change  COOE  (PLNARASCC) 
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SAMPLE  PROBLEM  2 
Swept  Wing  Configuration  With 
Transient  Heat  Conduction 


Sample  Problem  2  Input  Data 


0  0  0  0 

BRL  FLIGHT  CASE  (YUUA  TS-125  DEG-F.  T0-6P  OEG-F) 

TRANSIENT  COHOUCflON  SuLUTION  —  PLNRASCC  05  JANUAn  '  1984 

12.5  DEG  NOSE,  7  INCH  BOOT  <  BRL  Sample  planar  CASE>  Flnr  test 
01  Progrom  Conatonta  ond  Time  Inlormotion 
0.0  2.00 

4  0  15  0  1  e 


01 

02 


02 

03 


0  01 

0.01 

0  25 

0.25 

0.25 

2.00 

Env 

1  ronmen  t  Tob  1  e 

0  0 

1  0 

520. 

5259  . 

0  2 

1  0 

520. 

5184 

0  4 

1  0 

520 

5082 

0  6 

1  .0 

520 

4941 

0  8 

1  0 

520 

4793 

1  .0 

1 .0 

520. 

4636 

1  2 

1 . 0 

520 

4446 

1  .  4 

1  .0 

520 

4249 

1  .6 

1  .0 

520 

4035 

1  .8 

1  0 

520 

3839 

2.0 

1 .0 

520 

3629 

Sur 

f  oca  Geome  t  ry  and 

Grid  Si ze 

13 

0.2 

1  5 

0.4507 

e.e  1 

0.4627 

0  0684  1 

0.4975 

0  1286  1 

0.6074 

0.1953  1 

0.8 

.2380  1 

1  .0 

.2823  1 

1-249 

.3375  1 

1  55 

.3798  1 

19  .429  2 

50  .5836  2 

4.0  .7241  2 

5.5  0  9349  2 

70  1.1458  2 

34 

0.4507  0.0  1 

0.4627  0.0684  1 

0  4975  0  1286  1 

0.6074  0.1953  1 

0.8  .2380  1 

10  .2823  1 

1.249  .3375  1 

1  55  .3798  1 

1.875  .425  1 

1.94  .32  1 

1  76  .291  1 

1.7C  .226  1 

1.260  .1571  1 

1.268  .0775  1 

0  80  .0775  1 

0.80  0.0  1 


585 


75. 
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0  4507 

0.0 

1 

0.80 

0.0 

2 

0.80 

.0775 

2 

1  76 

.0775 

2 

1 .76 

.201 

2 

1 .94 

.32 

2 

1 .875 

.425 

2 

2.0 

.4431 

2 

4  0 

.7241 

2 

6.0 

1 .0052 

2 

7.0 

1 .1458 

2 

7.0 

0.0 

2 

0.60 

0.0 

2 

1 .268 

.0775 

3 

1  .268 

.1571 

3 

1  .76 

.2260 

3 

1  76 

.0775 

3 

1  .268 

.0775 

3 

03  7.0  0  0 

1 1 
.  1 

.  12 

60  22 


0  118© 

0  1 186 

0. 1186 

0  1186 

e  1  186 

.? 

1  186 

0  11 06 

0  1  1 86 

0.1186 

0  1  166 

0  1186 

ft 

1  186 

0.1186 

0.1186 

0.1186 

0.1186 

0  1  186 

C 

1  166 

0.1186 

0  1  186 

0.1186 

0  1  1B6 

0  1  186 

e 

1  186 

0  1  186 

0.1186 

0.1186 

0.1186 

0  1186 

e 

1  186 

0  1  186 

0.1186 

0.1186 

0  1186 

0.1186 

o 

1  1  96 

0  1 186 

0  1 186 

0  1 186 

0  1186 

0  1  186 

0 

1  186 

0.1186 

0.1186 

0.1186 

0  1 186 

e  1186 

0 

1  186 

0.1186 

0.1186 

0.1186 

0.1*86 

0  1  186 

c 

1  186 

0  1  186 

0  1186 

0.1186 

0  1  186 

0  1  186 

0  055 

0  055 

0.055 

0  055 

0  055 

0  056 

0  055 

0.055 

0.055 

0.055 

0.055 

0.055 

0  055 

0  055 

0.056 

0  055 

0  055 

0.055 

0.055 

0.055 

0  055 

0  055 

0  055 

0  055 

06  Material 

Prope  ''tie* 

Table 

1  e 

1 

10.0 

0  0 

418.08 

536  0 

0.0 

400  . 

.  063 

1  .  794E-02 

.  15 

500. 

.068 

1 .8056-02 

.15 

600. 

093 

1  7906-02 

.  15 

700. 

.098 

1 .758E-02 

.15 

800. 

.  103 

1 .7186-02 

.16 

900 

.  108 

1  670E-02 

.15 

1000. 

.112 

1  618E-02 

15 

1100. 

.117 

1  563E-02 

15 

-1  1200. 

.  122 

1  . 506E-02 

.  15 

2  0 

1 

10.0 

0  0 

490. 

536.0 

0.0 

46 


-1 


492. 

.  1  1 

7  360 E  -e' 

60 

672. 

.  1  1 

7  220E-02 

60 

1032. 

.  1  1 

6  940E-0? 

.  60 

1392. 

.  1  1 

6. 1  10E-05 

60 

3 

0 

1 

0.0 

0  0 

488  8 

536 

117.0 

3310 

50 

'28 

540 

.126 

7  56E-03 

5 

3250 

120 

7  S6E-0J 

5 

3320 

1932 

2  I0E-03 

5 

9000 

.  1932 

2  10E-0? 

5 

1 

2 

1  .  0E-05 

1 

3 

1  0E-04 

2 

3 

1 . 0E-06 

07  Surface  Pressure  Toble 

4  1 

60.  5 

+  1 

08  Shock  Shape  Table 

3 

09  Surface  Chemistry  Tobies 

1 

1  .0 


1  00 

00 

99  00 

662  210 

00 

1  00 

.00 

10  00 

659.444 

00 

1  00 

.00 

1  00 

656.666 

00 

1  00 

00 

0' 

653.888 

00 

1  00 

.00 

001 

600.000 

00 

1  00 

.00 

0001 

298.000 

00 

1.00 

.00 

00001 

200.000 

00 

2 

1 

0 

1  00 

00 

01 

2590.00 

.00 

1 .00 

.00 

001 

600.000 

00 

1.00 

.00 

0001 

298  000 

00 

1 .00 

.00 

00001 

200.000 

00 

3 

1  . 

0 

1  .00 

.00 

01 

2500.00 

.00 

1  .00 

.00 

001 

600.000 

.00 

1  .00 

.00 

0001 

298.000 

00 

1  .00 

.00 

00001 

200.000 

00 

99  556 

68  744 

1  A  r  •  •  . 

990+02 

89.850 

58  408 

i  A  :  '  •  i  • 

.  100+0? 

69  128 

58  087 

:  a*  i  :  . 

100+01 

80  405 

57  765 

:  a-;  •  ■ 

100-01 

7 4  555 

74.555 

1  A  , 

0. 

0 

\  a  r- 

4  2  130 

-42.130 

1  A  ! 

25  096 

525.096 

1  A  1 

7-»  555 

'’4.555 

1  A  1 

0 

0 

1  A 

4  2  130 

-  4  2  1  30 

1  A  1 

25  e96 

525  096 

1  A  j 

7  4  555 

’4.555 

1  A  1 

0 

0 

1  A  i 

42  130 

-42  .  130 

1  A  i 

-1 


-1 
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Sample  Problem  2  Output 

PLANAR  VERSION  BRl  IMPROVED  A8RES  SHAPE  CHANGE  COOE  (PLNARASCC) 


INPUT 


BRL  ELICHT  CASE  (YUMA  TS-125  DEG-E.  T0-60  DEG-E) 

TRANSIENT  CONDUCTION  SOLUTION  —  PLNRASCC  05  JANUARY  1984 

12  5  DEG  NOSE.  7  INCH  BOOY  <  BRL  SAMPLE  PLANAR  CASE>  PLNR  TEST 


-  GENERAL  PROGRAM  FLAGS  - 

(ENVIRONMENT  FLAG)  LC  -  4 

(SHAPE  CHANGE  FLAG)  ISS  -  0 

(OUTPUT  PRINT  FLAG)  IPRNT  -  1 

(TRANSITION  CRITERIA  FLAG)  NREYCR  -  5 

(BOOT  ANGLE  DEFN  FLAG)  IRON  -  0 

(CARBON  TRANS.  CRIT.  FLAG)  ICARB  -  1 

(NOSE  SHAPE  MODIFICATION  FLAG)  1UOO  -  0 

-  TIME  INCREMENT  INFORMATION  - 

INITIAL  TIME  (SEC)  0.0000  FINAL  TIME  (SEC)  2  0000 

OUTPUT  INTERVAL  -  0.0100  SEC  FROM  INITIAL  TIME  UNTIL  0  0100  SEC 

OUTPUT  INTERVAL  -  0  2500  SEC  FROM  0  0100  SEC  UNTIL  0  2500  SEC 

OUTPUT  INTERVAL  -  0.2500  SEC  FROM  0  2500  SEC  UNTIL  FINAL  TIME 


TIME  STEP  STABILITY  CRITERIA  IN  EFFECT 


MINIMUM  TIME  STEP  -  1 . 000E-06  SECONDS 


40 


planar  version  brl  improved  adres  SHAPE  CHANGE  COOE  (PLNARASCC) 


' 

H 

I  ^ 

-  CENERAL  ENVIRONMENT  - 


TIME 

PRESSURE 

[SEC) 

(ATM) 

e  eee 

1  .000 

0.200 

1  .000 

e.«00 

1  ,000 

0.600 

1  .000 

0.000 

1  .000 

1  .000 

1  .000 

1.200 

1  .000 

1 .400 

1  .000 

1 .600 

1  .000 

1  .808 

1  .000 

2.000 

1  .000 

TEMPERATURE 

VELOCITY 

(DEG  R) 

(EPS) 

520.00 

5259  00 

520  00 

5184.00 

520.00 

5082.00 

520.00 

4941  00 

520.00 

4793.00 

520  00 

4636  00 

520  00 

4446.00 

520.00 

4249.00 

520.00 

4035  00 

520.08 

3839  00 

520.00 

3623.00 

mllllll  "llutL’! Ill’Ll 


PLANAR  VERSION  BRL  IkPROVEO  ABRES  SHAPE  CHANGE  CODE  (PLNARASCC) 


PLANAR  VERSION  BRL  IMPROVED  ABRES  SHAPE  CHANGE  COOE  (PLNARASCC) 


- IMPLICIT  NOOE  SPACING - 

NOOE  THICKNESS  IN  INCHES 


NOOE  NO 

1  1 

2 

3 

4 

5 

6 

7 

8 

9  10  11 

BOOY  PT  NO. 
1 

i  0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.8120 

0.0120 

0.0120 

0  0120 

2 

1  8.0128 

0.0120 

0  0120 

0.0120 

0.0120 

0.0120 

0.0(20 

0  0120 

0  0120 

0.0120 

3 

1  0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.012e 

0.0120 

0.0120 

0.0120 

0.0120 

4 

1  0.0120 

0.0120 

0  0120 

0  0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.012O 

5 

1  0.0120 

0.0120 

0  0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0  0120 

6 

1  0.0120 

0.0120 

0  0120 

0.0120 

0  0120 

0.0120 

0.0120 

0  0120 

0.0120 

0.0120 

7 

1  0.0120 

0.0120 

0  0120 

0  0120 

0.0120 

0.0128 

0.0120 

0  0120 

0.0120 

0  0120 

8 

1  0  0120 

0.0120 

0  0120 

0  0120 

0  0120 

0.0120 

0.0120 

0  0120 

0  0120 

0  0120 

9 

1  0.0120 

0  0120 

0  0120 

0.0120 

0.0120 

0.0120 

0  0120 

0  0120 

0  0120 

0.6120 

ie 

1  0  0120 

0.0120 

0  0120 

0.0128 

0  0120 

0.0120 

0.0120 

0  0120 

0.0120 

6.0120 

n 

1  0.0120 

0.0120 

0  0120 

0.0120 

0.0120 

0.0120 

0.0120 

0  0120 

0  0120 

0  0120 

12 

1  0.0120 

0.0120 

8  0 1  20 

0.0120 

0.0120 

0.0120 

0.0(20 

0  0120 

O.012O 

0.9126 

13 

1  0.0120 

0.0120 

0  0120 

0  0120 

0.0120 

0.0120 

0  0120 

0  0120 

0.0120 

0  0120 
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-  EXPLICIT  GRID  GEOMETRY  - 

NLAfiER  OF  COLUMNS  -  60 
NUMBER  OF  ROWS  -  22 


VARIABLE  GRID  SPACING  WITH 
X  GRID  SPACING  - 


0.1 1660 

0.11860 

0.11860 

0.11860 

0.1 1860 

0.11660 

0.11860 

0. 11860 

0.1 1860 

0.  11860 

0.1 1860 

0.1 1860 

0.11860 

0.11860 

0.  1 1860 

0.1 1860 

0.11860 

0.11860 

0.11860 

0  11860 

0.11660 

0.11860 

0.11860 

0.11660 

0.11860 

0.1 1860 

6.11860 

0  11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0. 11860 

011 860 

0.11860 

0.1 1860 

0.  11860 

0.11 860 

0. 11860 

0.11860 

0.11860 

0.11860 

0  1 1860 

0.11860 

0.11860 

0.11860 

0.  11860 

0.11860 

0.1 1860 

0  1 1860 

0.11860 

0. 1 1860 

0  1 1860 

0.11860 

0.11860 

0.11860 

0.11860 

Y  GRID  SPACING  - 

0  05500 

0.05500 

0.05500 

0 . 05500 

0.05500 

0  05500 

0.05500 

0.05500 

0  05500 

0  05500 

0.05500 

0.05500 

0.05500 

0  05500 

0  05500 

0.05500 

0.05500 

0.05500 

0  05500 

0.05500 

0.05500 

INITIAL 

TEMPERATURE  OF  MOOEL  -  ! 

585.0  DEG  R 

MAXIMUM  DESIRED  SURFACE 

TEMPERATURE  RISE  BETWEEN 

TIME  STEPS  ' 

-  75.0  DEC  R 

MINIMUM  EXPLICIT  NODAL 

SPACING  USED 

IN  TIME  STEP 

COMPUTATION 

-  0.0550  INCH 
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- — MATERIAL  PROPERTIES - 


MATERIAL  NUMBER  1 


-  SURFACE  ROUGHNESS  - 

ROUGHNESS  HEIGHT  FOR  TRANSITION  K-LAM  .  US  000  (MIL) 

ROUGHNESS  HEIGHT  FOR  TURBULENT  HEATING  K-TURB  -  0  000  (MIL) 

LAMINAR  HEATING  AUGMENTATION  FLAG  JROUGH  -  1 


-  PARAMETERS  IN  BLOWING  CORRECTION  TO  TRANSFER  COEFFICIENTS  - 

LAMINAR  SHEAR  PARAMETER  (BLS)  -  0  5000 
LAMINAR  HEATING  PARAMETER  (BLH)  -  0  5000 
TURBULENT  SHEAR  PARAMETER  (RTS)  »  0  3500 
TURBULENT  HEATING  PARAMETER  (BTH)  -  0  3500 


-  THERMAL  PROPERTIES  - 

MATERIAL  DENSITY  (RHO)  -  418.08  (L0M/FT3) 

OATUM  TEMP  FOR  HEAT  OF  FORMATION  (TFO)  -  536  00  (DEG  R) 

HEAT  OF  FORMATION  (hFO)  -  0  00  (BTU/LBM) 


TEMPERATURE 

SPECIFIC  HEAT 

CONDUCTIVITY 

SENSIBLE 

EMI  Si  1 V] 

ENTHALPY 

(OEG  R) 

(BTU/LB-DEG) 

(8TU/FT-SEC-DEC) 

(BTU/LB) 

400.00 

0.0830 

0.0179400 

-1 1  81 

0  1500 

500  00 

0.8880 

0.0180500 

-3  26 

0  1500 

600 . 00 

0  0930 

0.0179000 

5.79 

0  1500 

700  00 

0 . 0980 

0.0175800 

15  34 

0  1500 

800  00 

0.1030 

0.0171800 

25  39 

0  1500 

900 . 00 

0  1080 

0.0167000 

35  94 

0  1500 

1000.00 

0  1  120 

0  0161800 

46  94 

0  1500 

1100.00 

0  1  170 

0  0156300 

58  39 

0  1500 

1 200 . 00 

0  1220 

0.0150600 

70.34 

0  1 500 

—  EROSION 

LAW  MATERIAL  FLAG 

_ 

NEROOE  -  0 
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MATERIAL  NUMBER  2 


-  SURFACE  ROUGHNESS  - 


ROUGHNESS  HEIGHT  FOR  TRANSITION 
ROUCHNESS  HEIGHT  FOR  TURBULENT  HEATING 


K-LAM  -  10  000  (MIL) 

K-TURB  -  0  000  (MIL) 


LAMINAR  HEATING  AUGMENTATION  FLAG 


J ROUGH  -  1 


-  PARAMETERS  IN  BLOWING  CORRECTION  TO  TRANSFER  COEFFIuiENTS  - 


LAMINAR  SHEAR  PARAMETER  (BLS) 
LAMINAR  HEATING  PARAMETER  (PLH) 
TURBULENT  SHEAR  PARAMETER  (BTS) 
TURBULENT  HEATING  PARAMETER  (BTH) 


0 . 5000 
0.5000 
0  3500 
0  3500 


-  THERMAL  PROPERTIES  - 


MATERIAL  DENSITY  (RHO) 
DATUM  TEMP  FOR  HEAT  OF  FORMATION  (TFO  j 
HEAT  OF  FORMATION  (HFO) 


490.00  ( LGM/FT  3  ) 
536  00  (DEG  R) 

0  00  (BTU/LBM) 


TEMPERATURE 

(OEG  R) 
492  00 
672  00 
1032  00 
1392  00 


SPECIFIC  HEAT 

(BTU/LB— OEG) 

0  1  100 
0.1100 
0.1100 
0  1100 


CONDUCTIVITY 

(BTU/FT-SEC-OFG) 
0  0073600 
0  0722060 
0  0694000 
0  0061100 


SENSIBLE 
ENTHALPY 
(BTU/LB) 
-4  84 
14  96 
54  56 
94  16 


EMissivm 


0  6000 
0  6000 
0  6000 
0  6000 


-  EROSION  LAW  MATERIAL  FLAG  - 


NEROOE 


0 
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♦••••••  MATERIAL  NUMBER  3  . . 


Planar  version  bri  improved  abres  shape  CHANCE  COOE  (PLNARASCC) 


CONTACT  RESISTANCES 


MAT  1 

MAT2 

RESISTANCE 
( F  t • • 2-S-0«gR/BTU) 

t 

2 

i .eeeeoe-es 

1 

3 

i . eeeeeE-64 

2 

3 

( . 00000 E-06 

PLANAR  SWEPT  WING  CASE 

swept  wing  ancle  -  so  00000  (oeg)  wedge 

ancle  - 

5  00000  (DEG) 

2-0  CVLINOER-WEDCE  PRESSURE  CORRELATIONS  ARE  USED 

MACH  NUflER  NORMAL  TO  LEADING  EDGE  -  2.34(82 
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— 

SHOCK  SHAPE  — 

— 

DIMENSIONLESS 

SHOCK  ANCLE 

OIMENSIONLESS 

Y-COORDINATE 

X-COOROINATE 

Y/RN 

(DECREES) 

X/RN 

0.600 

90.000 

0.000 

0.100 

88.772 

0.000 

0.200 

87.276 

0.000 

0  300 

85.819 

0.000 

0.400 

84.400 

0  000 

0.500 

83.018 

O.O0O 

0.600 

81 .673 

0.000 

0.700 

80  363 

0.000 

0  800 

79 .090 

0.000 

0.900 

77,850 

0  000 

1  .000 

76  645 

0  000 

1.100 

75.474 

0.000 

1  .200 

74.335 

0.000 

1  .  .100 

73.228 

0.000 

1 .400 

72.153 

0  000 

1  500 

71.109 

0  000 

1  600 

70.096 

0  000 

1  .700 

69.112 

0  000 

1  .800 

60. 157 

0  000 

1  .900 

67  230 

0.000 

2  000 

66,332 

0.000 

2.100 

65  460 

0  000 

2  200 

64.615 

0  000 

2.300 

03.797 

0.000 

2  400 

63  O03 

0  000 

2  500 

62.235 

0.000 

2.600 

61  .491 

0.000 

2  700 

60.770 

0  000 

2  800 

60.073 

0  000 

2.900 

59 , 397 

0.000 

3  000 

58  744 

0  000 

3.100 

58  11 3 

0  000 

3.200 

57.562 

0  000 

3  300 

56.911 

e  eeo 

3.400 

36.340 

0  000 

3.500 

55.78B 

0  000 

3  600 

55.255 

e  eoe 

3.700 

54.739 

0.000 

3  800 

54.242 

0  000 

3  900 

53.761 

0.000 

4.000 

53.297 

0.000 

4.100 

52.849 

0.000 

4.200 

52.416 

0.000 

4.300 

51 .999 

0.000 
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4.400 

51.596 

0.000 

4.560 

51.207 

0.000 

4.666 

50.832 

0.000 

4.766 

50.470 

0.000 

4.880 

50.120 

0.000 

4.900 

49.783 

0.000 

5.000 

49 . 458 

0.000 

5.100 

49  144 

0.000 

5.200 

48.841 

0.000 

5.300 

48.549 

0.000 

5.400 

48.267 

Ok  OAA 
.  Ovv 

5.500 

47.995 

0.000 

5.600 

47.732 

A  AAA 

D  .  VOW 

5.700 

47.478 

0.000 

5.800 

47.233 

0.000 

5.900 

46.996 

0.000 

6.000 

46.767 

0.000 

6.100 

46.546 

0.000 

6.200 

46 . 332 

e.eee 

6.300 

46.125 

0.000 

6.400 

45  925 

0.000 

6.500 

45.731 

0.000 

6.600 

45.543 

0.000 

6.700 

45.361 

6  000 

6.800 

45  184 

0.000 

6.900 

45.012 

0  000 

7.000 

44.846 

0.000 

7.100 

44 . 684 

0.000 

7.200 

44.527 

0  000 

7.300 

44.373 

0.000 

7.400 

44  224 

0.000 

7.500 

44.079 

A  AA" 

W  .  WWW 

7.600 

43  937 

O  ,  Wo 

7.700 

43 . 798 

0  000 

7.808 

43 . 663 

0  000 

7.900 

43.531 

0.000 

0.000 

43.401 

0  000 

8  100 

43.274 

0  000 

8.230 

43.113 

0.000 

8.399 

42.909 

0.000 

8  619 

42.653 

0  000 

8.904 

42.332 

0.000 

9.276 

41 ,933 

0  000 

9.758 

41.438 

0  000 

10.386 

40.827 

0.000 

11.201 

40 . 080 

0  000 

12  262 

39.198 

0.000 

13.641 

38.263 

0  000 

15.433 

37.455 

0  000 

17.763 

37.112 

0.000 

20.791 

36.538 

0 . 000 

24.729 

35.711 

0.000 

29.847 

34.637 

0.000 

36.501 

33.240 

0.000 

45.152 

31 .424 

0  000 

56  397 

30.406 

0  000 

71 .017 

30.406 

e  eee 
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- SURFACE  EQUILIBRIUM  DATA - 

MAT  -  1 

CMH  -  1.00000 

MASS  TRANSFER  COEF.  -  0.0000  CBM/FT.. 2-SEC  PRESSURE  -  1.0000  ATM 


T£W> 

360 . 0000 
536 . *000 
I860  0000 
1176  996* 
1181 .9988 
1186.9992 
1191 .9760 

BPR1M 

0  0000 
6.0001 

0  0010 

0  0100 

1 ,0000 
10.0000 
99.0000 

HCH 

-15.2280 

0  0362 

56  1020 

67  5933 

6B  <909 

68  788* 

69 . 383* 

TSEN 
-75  83*0 
0.0000 
13*. 1990 
1039770 
10*. 5566 
106. 13** 
105.7392 

TCHEM 

75  83*6 
0.0000 
-13*  2771 
-159.4928 
-196  7961 
-525.1900 
-3778. *279 

SPECIE  i 

AIR 

AIR 

AIR 

AC*  1 

AC*  1 

AC*  1 

AC*  1 

MAT 

CMH 

-  2 

.  1.00000 

MASS  TRANSFER 

COEF.  «  0  0000 

LBM/FT«»2- 

-SEC 

PRESSURE  - 

1  0000  ATM 

TEMP 
360.0000 
536 . *000 
1080.0000 
*5000000 

8PR1M 

0  0000 

0  0001 
0.0010 

0  0100 

HCH 

-19  3600 
0.0*40 

59 . 8*00 
*36.0*00 

TSEN 

-75.83*0 

0  0000 
13*. 1990 
9*5. 1728 

TCHEM 

75  6346 

0  0000 
-13*  273* 
-95026*1 

SPECIE 

AIR 

AIR 

AIR 

AIR 

MAT  -  3 

CMH  -  1 . 00000 
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NOMENCLATURE 


H  convective  heat  transfer  coefficient  in  Figures  2-5  through  2-11 

A  swept  wing  angle  measured  from  a  normal  to  the  flow  direction 

Me  frees^ream  Mach  number 

P  static  pres  ,re 

P0  stagnation  pressure 

Re  Reynolds  number 

R-j  nose  radius  on  cylinder/wedge  swept  wing 

S  distance  along  surface  measured  normal  to  the  wing  leading  edge 

T  static  temperature 

T0  stagnation  temperature 

X  axial  distance  measured  nor  '  to  the  wing  leading  edge 

ew  aft  wedge  angle  on  cylinder/wedge  swept  wing 
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